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O . ABSTRACT 
(N 

Weyl's idea on scale invariance was resurrected by Cheng in 1988. The requirement of local scale 
invariance leads to a completely new vector field, which we call the "Cheng- Weyl vector field" . The 
Cheng- Weyl vector field couples only to a scalar field and the gravitational field naturally. It does 
not interact with other known matters in the standard model of particle physics. In the present 
00 ' work, the (generalized) Cheng- Weyl vector field coupled with the scalar field and its cosmological 

application are investigated. A mixture of the scalar field and a so-called "cosmic triad" of three 
mutually orthogonal Cheng- Weyl vector fields is regarded as the dark energy in the universe. The 
cosmological evolution of this "mixed" dark energy model is studied. We find that the effective 
equation-of-state parameter of the dark energy can cross the phantom divide Wde = — 1 in some 
cases; the first and second cosmological coincidence problems can be alleviated at the same time in 
this model. 
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I. INTRODUCTION 

Since the discovery of present accelerated expansion of our universe p], 0, H, 0, [E @, 0] , dark energy Q 
has been an active field in modern cosmology One of the puzzles of the dark energy problem is the 
(first) cosmological coincidence problem, namely why does our universe begin the accelerated expansion 
recently? why are we living in an epoch in which the energy densities of dark energy and dust matter are 
comparable? In order to give a reasonable interpretation to the (first) cosmological coincidence problem, 
many dynamical dark energy models have been proposed, such as quintessence U 0] ■> phantom (ill . Il2l . Il3j , 
k-essence [HI, EH etc. 

The equation-of-state parameter (EoS) of dark energy Wde = Pde/ Pde plays a central role in observational 
cosmology, where pde and pde are its pressure and energy density respectively. Recently, by fit ting the 
observational data, marginal evidence for Wd e (z) < —1 at redshift z < 0.2 ~ 0.3 has been found [l6l.ll7l [l8j. 
In addition, many best-fits of the present value of Wde are less than — 1 in various data fittings with different 
parameterizations (see [l9[ for a recent review). The present observational data seem to slightly favor an 
evolving dark energy with Wde crossing —1 from above to below in the near past [l7| . Obviously, the EoS 
of dark energy Wde cannot cross the so-called phantom divide Wde = — 1 for quintessence or phantom alone. 
Although at first glance, it seems possible for some variants of k-essence to give a promising solution to 
cross the phantom divide, a no-go theorem, shown in [2(il |. shatters this kind of hope. In fact, it is not 
a trivial task to build dark ener gy m odel whose EoS can cross the phantom divide. To this end, a lot 
of efforts 0000 01S, 0000 [fii] have been made. 
To name a few, quintom model, string theory inspired models, vector field models, crossing the phantom 
divide in the braneworld models, scalar-tensor models 0], etc. However, to our knowledge, many of those 
models only provide the possibility that Wde can cross —1. They do not answer another question, i.e., why 
does crossing the phantom divide occur recently? why are we living in an epoch Wde < ~1? This can be 
regarded as the second cosmological coincidence problem 0, 0] . 

Although cosmological observations hint that Wd e < —1, however, there is a subtle tension between 
observations and theory. For the canonical scalar phantom model , the universe has an inevitable fate 
of big rip [l2|], and the instability is inherent [lj^]. For the k-essence model with EoS less than —1, the 
spatial instabilities inevitably arise too 0] (see also [HI)- In a more general case, it is argued that there 
is a direct connection between instability and the violation of null energy condition (NEC) 0,0]. Some 
models can evade this result, at the price of the lack of isotropy of the background and the presence of 
superluminal modes 0, 0] • Recently, two seemingly viable models that violate NEC without instability 
or other pathological features have been proposed in different ways 0, 0]. In particular, in [46], a scalar 
field coupled with a vector field is used; and the effective Lagrangian explicitly depends on the vector 
filed Afj,, which avoids one of the assumptions of 0, 0] that the effective Lagrangian only depends on 
F^ u = d^A u — d„A^ and has no dependence on the vector filed A^ itself. 

Motivated by the work of [4a], it is interesting to study the case of a scalar field coupled with a vector 
field. In fact, the vector field has been used in modern cosmology in many cases, see e.g. 0, 0| and 
references therein. It is worth noting that comparing to the ones investigated in e.g. 0,0], the vector 
fields used in 0, S3] have fairly different motivations and forms. In 0], a single dynamical scalar field 
is coupled with an a priori non-dynamical background vector field with a constant zeroth-component. In 
that model, the effective EoS can cross the phantom divide w = — 1. However, the appearance of the a 
priori non-dynamical vector field has no clear physical motivation. 

In 0, 0], Weyl's old idea of scale invariance 0, 0] was resurrected by Cheng in 1988, almost 60 
years later (See [52j for an independent rediscovery). The requirement of local scale invariance leads to 
the existence of a completely new vector field, which we call the "Cheng- Weyl vector field" throughout 
this paper, in honor of the proposer Cheng and Weyl (a great mathematician and physicist [53]). The 
Cheng- Weyl vector field only couples to the scalar field and the gravitational field. It does not interact 
with other known matters in the standard model of particle physics, such as quarks, leptons, gauge mesons, 
and so on. In particular, it has no interaction with photons and electrons. So, it is "dark" in this sense. 
As mentioned above, the fact that Cheng- Weyl vector field naturally couples to the scalar field makes it 
very interesting, especially when the scalar field is considered as dark energy candidate. Required by the 
local scale invariance, the potential term of the scalar field has to be of </> 4 form, while the coupling form 
between the Cheng- Weyl vector field and the scalar field is fixed also, and the form is different from the 
ones of . 0] . Interestingly, the effective Lagrangian also explicitly depends on the Cheng- Weyl vector 
field itself naturally. In Sec. [Til we will give a brief review of the work of Cheng 0, 0], in which the 
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Cheng- Weyl vector field was proposed. 

In the present work, the (generalized) Cheng- Weyl vector field coupled with a scalar field and its cos- 
mological application are investigated. We regard a mixture of the scalar field and a so-called "cosmic 
triad" of three mutually orthogonal Cheng- Weyl vector fields as the dark energy in the universe. We derive 
the effective energy density and pressure of the "mixed" dark energy, and the equations of motion for the 
scalar field and the Cheng- Weyl vector field respectively. The cosmological evolution of this "mixed" dark 
energy is studied. We find that the effective EoS of dark energy can cross the phantom divide Wd e = — 1 
in some cases; the first and second cosmological coincidence problems can be alleviated at the same time 
in this model. 

This paper is organized as follows. In Sec. [Til we wm briefly present the main points of the Cheng- Weyl 
vector field proposed in [48|, [4t|. In Sec. IIII1 the effective energy density, pressure, and the equations of 
motion are obtained. In Sec. IIVI The cosmological evolution of the "mixed" dark energy is investigated 
by means of dynamical system |60l |; the first and second cosmological coincidence problems are discussed. 
Finally, a brief conclusion is given in Sec. [Vl 

Throughout this paper, we use the units % = c = 1 and the notation k 2 = 8irG, and adopt the metric 
convention as (+,—,—,—). 



II. 



THE CHENG- WEYL VECTOR FIELD 



Following [48|, |49( , here we give a brief review of the so-called Cheng- Weyl vector field. The arguments 
are based on the local scale invariance. It is important to distinguish the scale invariance from the gauge 
invariance. The scale invariance is the invariance of the action under the change of the magnitude rather 
than the phase of the fields. To be definite, let us consider the distance between two neighboring spacctime 
points, ds 2 = g^ v dx^dx v . We change the scale of the distance, for instance, changing the unit of length 
from meter to inch. With this change, the distance remains the same, but is measured in a different 
unit That is, 



9pv * 9pv 



(1) 



where A is a constant for the global scale invariance, and is a function of space and time for the local 



scale invariance. Then we have ds 2 — > ds 2 = g liv dx ii dx v = A 2 ds 2 , and g 



9 



A 



= A 4 ^/[g|, where g is the determinant of the metric <? M „. So, the action S = J d 4 
invariant under the scale transformations, provided that the Lagrangian density satisfies 



And, 
\g\C is 



C^C = A~ 4 C. 



(2) 



In this case, the forms of all equations in the theory remain the same. 

Let us first see the case of the global scale invariance, i.e. A is a constant. The Lagrangian density of a 
scalar field is given by 



1 



g^d^du^ - A0 4 



(3) 



where A is a dimensionless constant. The Lagrangian density of a gauge meson is 



^y y 1 pm 



(4) 



where F M „ = d^Ay — <3„A M , or — d^W" — d v W® — gf abc W^W£ (here g is a coupling constant) for the 
Yang- Mills theory. The Lagrangian density for a fermion \? coupled with the electromagnetic field and the 
gravitational field is 



bv 



(5) 



9 pa {dp9a V + dv9op ~ d a g^)/2, a ab = ( 7 Q 7 b - 7 V)/4, and e« is the tetrad satisfying 



where 

Vab^st = 9pv It is easy to verify that the above Lagrangian densities satisfy Eq. ^ under the scale 
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transformation Eq. ([T]) and 

«^ = A"V, A li ^A lt = A li , WZ^WZ = W£, eJ-e^AsJ, * _> $ = A" 3 / 2 *. (6) 

Next we consider the case with the theory being scale invariant locally, i.e. A is a function of space and 
time. Similar to the well-known arguments used to deduce the existence of gauge fields, one can find that 
a completely new vector field ±I M , namely the so-called Cheng- Weyl vector field, is required by the local 
scale invariance, while the replacements 

d^gu P -> (fy + 2/H M ) g vp , d^g"" -> (0„ - 2/ILJ 

^-(^-/n^)^, fy* -> (fy - §/n„) *, (7) 

are also required in these Lagrangian densities, where / is a dimensionless constant. One can verify that 
these Lagrangian densities with the replacements Eq. satisfy Eq. |2j under the scale transformation 
Eqs. (HD, ©and 

n^-ffi^ = n /1 --a |t inA. (8) 

With the replacements Eq. Q , the Lagrangian density of the scalar field Eq. ((3| becomes 

^gr^-fnjtifr-fTLJt-Xif. (9) 

Thus, the scalar field is coupled with II M naturally by the scale invariance. However, with the replacements 
Eq. 0, the Lagrangian densities of the gauge meson and the fermion, i.e. Eqs. (H)) and (0 respectively, 
need not to be altered, since the terms involving IT M completely cancel one another [H, |4§|. Therefore, 
we conclude that the gauge meson and the fermion do not couple with ILj. With identical arguments, 
the quarks and leptons etc. do not couple with 11^ as well. After all, we would like to mention that the 
Lagrangian density of 11^ itself [48| , 

- ^V ff W, (io) 

also satisfies Eq. under the transformations Eqs. |T]), © and (©, where 

y M „ = 9 M n„ - d v n^. (ii) 

We close this section with some remarks. First, required by the local scale invariance, the potential term 
of the scalar field has to be of <fi 4 form, as in Eqs. ^ and @. Second, it is easy to see from Eq. (O that 
the coupling form between the Cheng- Weyl vector field 11^ and the scalar field 4> is fixed naturally. Note 
that the form is quite different from the ones considered in [13, l46j|. For more details on the Cheng- Weyl 
vector field, please see (4|| [4j| . In addition, one may also refer to |54l. [55j for relevant papers on the scale 
invariance. 



III. APPLYING THE CHENG- WEYL VECTOR FIELD TO COSMOLOGY 



In modern cosmology, the scalar field is used extensively. Actually, the scalar field is one of the leading 
dark energy candidates. If the nature respects the local scale invariance, the so-called Cheng- Weyl vector 
field must exist, and couples to the scalar field inherently. If the scalar field is indeed the cause driving the 
accelerated expansion of the universe, we argue that the dark energy should be a mixture of the scalar field 
and the Cheng- Weyl vector field, which can be considered as the partner of the scalar field. This seems 
quite plausible when the fact that the Cheng- Weyl vector field does not interact with other known matters 
(so, it is "dark" to them) is taken into account. Therefore it is quite interesting to study the cosmological 
consequence of the Cheng- Weyl vector. 
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We begin with the action 



S = S. 



grav 



■S r 



where S grav and S m are the actions for gravitational field and matters respectively, and 



5V.w — / d 



(12) 
(13) 



Naively, one may write the Lagrangian density £ cw as 

1 



^ff" (d^ - jTL» {d v cj> - fU v <f,) - A0 4 - -g^g^Y^Y^, 



directly from Eqs. {9|) and (|10|) . In order to be compatible with homogeneity and isotropy, the IT^ can be 
chosen as IT M = (IIo, 0, 0, 0), where IIo = Ho(t) only depends on the cosmic time t. However, in this case 
= 0. From this Lagrangian density £ cw , for the case of homogeneous (j>, one finds that IIo = 4> — 
which is not dynamical, and the Lagrangian density £ cw is zero actually. Thus, unfortunately, the naive 
approach is not viable. 

Enlightened by the work of [47[ (see also [HI), we can describe the dark energy as a mixture of a scalar 
field and a so-called "cosmic triad" (in the terminology of [47) ) of three mutually orthogonal Cheng- Weyl 
vector fields. In this case, the Lagrangian density £ cw is given by 



3 

E 

a=l 
3 



■5 



a=l 



/ j L cw' 



4 - A0 4 - efg^WJd^ + -g^fWJLW 2 - ^g^Y^Y^ 



(14) 



where 



Y a = 8 TT a - 8 TT a 

1 n,u — u n LL i> u v s - s -ui 



(15) 



e, T], A and / being dimensionless constants. Latin indices label the different Cheng- Weyl vector fields 
(a, 6, . . . = 1, 2, 3) and Greek indices label different spacetime components (/x, v, , . . = 0, 1, 2, 3). Actually, 
the number of Cheng- Weyl vector fields is dictated by the number of spatial dimensions and the requirement 
of isotropy [H, The Latin indices are raised and lowered with the flat "metric" 5 a b- It is worth noting 
that the Lagrangian density £ cw in Eq. p4|) satisfies the requirement of the local scale invariance. Note 
that in Eq. (fT4|) . we have generalized the original scalar field to include the cases of quintessence (e = +1) 
and phantom (e = —1), while the Lagrangian density for the Cheng- Weyl vector fields has also been 
generalized by introducing the constant r\. 

Varying the action (fT3"|) with Eq. (fB|) . one can get the energy- momentum tensor of the "mixed" dark 
energy as 



where 



(cw)t 



« = e ( ~fl 



_2_^ 



E 



-9m 



/•(<*) 



n cw 
SgM" 



1 



(16) 



(17) 



From the action (TT3"|) with Eq. (O, one can also obtain the equations of motion for <f> and II" namely 



-.9^ (-4A0 3 - efg^d^ + efg^U^ 

a=l 



(18) 



G 



and 



vd» W~gg^ a 9 vfi y^) = tV~9 {fsT^ ~ <T/ 2 n> 2 ) , (19) 

respectively. We consider a spatially flat Friedmann- Robertson- Walker (FRW) universe with metric 

ds 2 = dt 2 - a 2 (t)dx 2 , (20) 

where a(t) is the scale factor. In this work, we assume the scalar field is homogeneous, namely cf> = (f>(t). 
Similar to [13], an ansatz for the Cheng- Weyl vectors is 



n b = s b u n(t) ■ a(t). 



(21) 



Thus, the three Cheng- Weyl vectors point in mutually orthogonal spatial directions, and share the same 
time-dependent length, i.e. IF 2 = -g^Nfil® = FI 2 (i). Hence, the equations of motion JTBJ and ([15)1 
become, respectively, 



U + 3H<pj + e/ 2 n 2 + 4A0 3 = 0, 



and 



V 



II + 3HU + H 



n 



e / 2 n> 2 = o, 



(22) 



(23) 



where H = a/ a is the Hubble parameter, and a dot denotes the derivative with respect to the cosmic 
time t. From Eqs. (H]), dHJ) and $T7J), we find that 



Pew 



= (cw) T = 3 ^2 + 3A04 + 3 ef 2 u 2^2 



2 " 2 
where p CVJ is the energy density of dark energy, while 

(cw) T = Q) 

and 



3 / . 2 



(24) 
(25) 



(cw) 



rpi \ A {a)rpi 

3 / j J 



(26) 



where 



(a) rpi 



l^_A0 4 -^/W + |(ii + i/n)' 



ef 2 U 2 (j} 2 -r)(tl + HIiy 



6\S a v (27) 



It is worth noting here that in fact, even adopting the "cosmic triad" of three mutually orthogonal Cheng- 
Weyl vector fields, Lagrangian (fT4"|) is not invariant under 50(3) rotation in the internal space. Thus, the 
energy-momentum tensor is not strictly diagonal. One can see this point by noting that ' cw 'Tq =/= 0, 
due to the second term in the right hand side of Eq. (|17|). To overcome this inconsistence with the 
isotropy, the spatial volume-average procedure has to be employed here as done in [56l H3] • In those two 
papers the authors considered the nonlinear electromagnetic field as the source driving the accelerated 
expansion of the universe. There, in order to obtain an energy-momentum tensor consistent with the 
FRW metric, the spatial volume-average procedure in the large scale [H, l59l ] has been used. By using this 
procedure, the spatial volume-averaged non-diagonal components of the energy-momentum tensor become 
zero, namely ^ CW ^T J ) = 0, while the diagonal components are kept unchanged. Therefore, in our model, 
the energy-momentum tensor is a spatial volume-averaged one on the cosmological scale. In this way, the 
energy-momentum tensor is compatible with isotropy. 

The corresponding pressure of dark energy is given by 



(cw) T *. = _ £( 



3A0 4 - -f 2 n 2 4> 2 



HU 



(28) 
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The Friedmann equation and Raychaudhuri equation read, respectively, 



H 



-Ptot 



(Pcw + Pm) ■ 



(29) 



and 



n 2 n 2 

H = (ptot +Ptot) = ^ ^cv 



Pn 



Pc 



Pm) , 



(30) 



where p m and p m are the pressure and energy density of the matters, respectively. 
From Eqs. {53]) and pg|l. we obtain 



Pew + Pew = 3e 2 + e/ 2 nV + 2?7 (rr + HU 



(31) 



Obviously, the EoS of dark energy u> cw = P C w/Pcw i s always larger than —1 for the case of e > and 
rj > 0, while u> cw < —1 for the case of e < and n < 0. Crossing the phantom divide is impossible for 
both cases. However, for the case of e and rj having opposite signs, u> cw can be larger than or smaller than 
— 1. Of course, for this case, crossing the phantom divide is possible. 
In addition, from Eqs. (f24|) and (|28|) . we have 



Pew + 3H (Pc 



30 



e (d> + ZH<j> 

(ri + Huj { 



e/ 2 n 2 + 



V 



n + 3i?n + 2H 2 + h)u 



(32) 



Noting that a/a = H 2 + H and the equations of motion for <j> and II, i.e. Eqs. (|22p and (|23|) . it is easy to 
see that the energy conservation equation of dark energy holds, namely, p cw + 3H (p cw + P cw ) = 0. 



IV. DYNAMICAL SYSTEM AND COSMOLOGICAL EVOLUTION 



In this section, we investigate the cosmological evolution of the "mixed" dark energy by means of 
dynamical system [60j. Our main aim is to see whether this model can alleviate the coincidence problems. 
Like many considerations in the literature, we allow the existence of interaction between the dark energy 
and the background matter (usually the cold dark matter). The cases of the scalar field and vector field 
interacting with background matter are studied extensively, see, for examples, [36, 61, 62, 63, 64, 65, 66H67l|. 
Although the Cheng- Weyl vector field does not interact with the known matters in the particle physics 
standard model, nothing precludes the possibility of the Cheng- Weyl vector field interacting with the cold 
dark matter, since the nature of cold dark matter is also unknown so far. 



A. Dynamical system 

We assume that the dark energy and background matter interact through interaction terms C and Q, 
namely 

Pcw +3#(p cw +P CW ) = -C-Q, (33) 
p m + 3H (p m +p m )=C + Q, (34) 

which keep the total energy conservation equation p to t + 3 if (ptot + Ptot) = 0. The background matter is 
described by a perfect fluid with barotropic equation of state 

Pm = WmPm = (j ~ l)pm, (35) 

where the barotropic index 7 is a dimensionless constant and satisfies < 7 < 2. In particular, 7=1 
and 4/3 correspond to dust matter and radiation, respectively. Due to the interaction terms C and Q, the 
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equations of motion for </> and II, namely Eqs. (|22|) and ([23| . are altered. Seeing from Eqs. (|32| and ([33]) . 
they are given by 



30 



n + nuj { 



(d> + 3H<p \ + efH?<f> + 4A</> 3 



-c, 



= -q, 



n + 3i?n + [2H Z + Hju 

respectively. Following [6fJ, |61|, |62j, |63|, [U, [65| , we introduce following dimensionless variables 

6 kU 



K(p K,H 

y ■ 



Ky'Pm (f> 

z = — -= — , u = — , V = — =. 
V3H H' V2 



(36) 



(37) 



With the help of Eqs 
a dynamical system }60j . i.e 



V2H 1 y "V2i?' 
((301) i C3) an( i CHI) , the evolution equations (|3l)|) and (|3"4"|) can be rewritten as 



(6i -3)x- 2V2 Xe- 1 uQ 2 - V2 f 2 v 2 u 3 ©3 1 - C u 



y' = (0! - 3) y + (0! - e ?T7 V - 2) « - Qi, 



z = 

u = u 
v' = y, 



0i - -1JZ + C2 + Q2, 
0! +V2xuQ^\ , 



where 



and 



C x = 



K 2 C 

6eH 3 a 



0i 



02 



IP 



n 2 Q 



6i]H 3 (y + v) ' 



r _ zC _ zQ 



2Hp, 



2Hp„ 



3ex 2 + ef 2 u 2 v 2 + 2rj(y + v) 2 + -73 
2, ,2, ,2 1 „^2 i V 2 



1 — z 2 — ex — ef 2 u 2 v — rj(y + v) 
A 



a prime denotes derivative with respect to the so-called e-folding time iV = In a. 

The fractional energy densities of the background matter and dark energy are given by 



2 

K Pm 2 

- Z , 



3H 2 



and 



a 



K Pc 



ex 2 + Xk 2 H 2 u 4 + ef 2 u 2 v 2 + r}(y + v) J 



cw- 3H 2 

respectively. On the other hand, from Eq. l[2"9"]). one has 

Hence, from Eqs. (|47|) and (|48| . one can find out 

kH = — , 

where ©2 is given byEq. (|45[l . The effective EoS of the whole system is 

Ptot 



"eff 



Ptot 



(38) 
(39) 

(40) 

(41) 
(42) 



(43) 

(44) 
(45) 

(46) 

(47) 
(48) 
(49) 
(50) 



where w cw s 
respectively. 



Pew/ Pew ano - w m = Pm/ Pm = 7 — 1 are the EoS of dark energy and background matter, 
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B. Interaction terms and critical points 



In this subsection, we obtain all critical points of the dynamical system (|3"5)) -([4"2" j) . A critical point 
(x,y, z, u, V) satisfies the conditions x' — y' — z' — u' = v' — 0. Before giving the particular interaction 
terms C and Q, let us first find the general features of the critical points of dynamical system (f3"5|l -([4"2" |) . 
From Eq. (|42l) and v' = 0, it is easy to see that 

V = 0. (51) 

If this dynamical system has some critical points, their corresponding x, y, z, u, and v should be constants. 
Therefore, from Eqs. (|49|) and (|45l) . the corresponding Hubble parameter H = H = const. From Eq. |44|) . 

9i = (52) 
follows. Substituting into Eq. f4"T]) . u' — requires 

X = 0, (53) 
since u 2 ^ 1 = (nH)- 1 ^ 0. Hence, from Eqs. (JSU), gl), jEE]) and $53$, we have 

(efu 2 + 2rj) v 2 + ^jz 2 = 0. (54) 

So, ef 2 u 2 + 2r/ < is required for non- vanishing real z and v. By using Eqs. (|38p -([40 |) and Eqs. ([5T|) - ((54|) . 
x' = y' = z ' = become, respectively, 



V2u(e9 2 ) 1 (2-jz 2 )+Ci = 0, (55) 



ef 2 u 2 + 2n)v + 7781 = 0, (56) 
C2 + Q2- \iz = 0, (57) 



where 



e 2 



, 1/2 

1 + (§7 -l)z 2 + nv 2 
A 



(58) 



which comes from Eqs. (j4"5j) . (l5T|) . ([53^) . and ([54^1 . Then, one can find out the remaining z, u and v from 
Eqs. ([54|) -([57 |) . Obviously, only three of them are independent of each other. 

So far, the above results are independent of particular interaction terms C and Q. To find out z, u and v, 
we have to choose proper C and Q here. The interaction forms extensively considered in the literature 
(see [U S, H| lellfil lei, M, M, m, [63 for instance) are 

C cx Hp m , Hp tot , Hp 

cw 1 K Pm ( f l ! ■ ■ ■ 

Q(xHp m , Hptot, Hp cw , Kp m n, ... 

Noting that Eqs. ([53")) . (|3T|) and the definition of C\ in Eq. JM]), we have to choose 

C = anp m 4>, (59) 

to avoid the divergence of C\ in Eq. (|55|) . where a is a dimcnsionless constant. In this case, from Eq. (|43|) . 
one has 

From Eqs. (O and (HO]), we find that C 2 = in Eq. ((571). Noting that Eqs. ![5T|). (|53l> . (gU) and the definition 
of Q2 in Eq. (f4"3"|) . we cannot choose Q cx i//9 m or Kp m Il, to avoid the solution z = from Eq. (|5T|) . since 
our main aim is to alleviate the coincidence problems. Therefore, we choose Case (I) Q = 3(3Hp cw or 
Case (II) Q = SaHp to t, where (3 and a are dimensionless constants. 
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Case (I) Q = 3f3H Pcw 

In this case, from Eq. (|43[) . one has 



Noting that C 2 = in Eq. (57]), we find out 



/3 + 7 



(62) 



One can check that Eq. (|56|) is equivalent to Eq. (|57|1 for this case. Then, one can find out u and v from 
Eqs. ([54|) and ((55|) . by using Eqs. (fSTjl . (|53|) and (|62|) . We do not present them here, since the final results 
are involved and tedious. One can work them out with the help of Mathematica. Instead we would like to 
give several particular examples to support our statement. Example (1.1), for parameters 7 = 1, e = 1, 
r) = -1, a = 5, /3 = 3/5, A = 0.1 and / = 5, we have u = -0.246341, v = ±1.07927. Example (1.2), for 
parameters 7 = 1, e = 1, 77 = —1, a = 5, p = 1/2, A = 0.1 and / = 5, we find that u = —0.249803, 
v = ±1.06605. 

From Eqs. (gBJ), (03), and ([52}. the fractional energy densities of background matter and dark energy 
are given by 

ft- = pT~y "cw = J^, (63) 

respectively. For reasonable f2 m and r2 cw , it is easy to see that P > is required. As mentioned above, at 
the critical point (x, y, z,u,v), the Hubble parameter H = H = const. From Eq. (|30[) . this means 

w eff = - 1 - ( 64 ) 
From Eq. we find that the EoS of the dark energy is given by 

w cw =-l-/3. (65) 

Obviously, w cw < — 1. 



Case (II) Q = 3aHp to t 



In this case, the corresponding Qi and Q2 read 



^ 2?7(y + u)' ^ z 2 
respectively. Solving Eq. (57J with G% — 0, we get 



3ff -2 = Jw- 1 , (66) 



(67) 



Again, one can check that Eq. (|56j) is equivalent to Eq. (|57j) for this case. Then, one can find out u and 
S from Eqs. (|5i|) and (55]) , by using Eqs. (|5T|) . ([B^|) and (|6T|) . Once again, we do not present the long 
and involved expressions here. We only give some particular examples. Example (II. 1), for parameters 
7 = 1, e = 1, 77 = -1, a = 4, a = 1/3, A = 0.1 and / = 3, we find that u = -0.410797, v = ±1.01934. 
Example (II. 2), for parameters 7 = 1, e = 1, 77 = —1, a — 5, a — 0.3, A = 0.1 and / = 7, we have 
u = -0.180804, v = ±1.06307. 
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From Eqs. (|46| . ([48]) . and (|67| . the fractional energy densities of background matter and dark energy 
are given by, respectively, 

fim = -, f! C w = 1 --. (68) 

7 7 

which requires < a < 7. Following a similar argument, we have w e ff = — 1 also. And then, from 
Eq. ([5"Uf , we find that the EoS of the dark energy is given by 

t" ow =-l — , (69) 

7 — a 

which is also smaller than —1. 



C. Stability analysis 



In this subsection, we discuss the stability of these critical points. An attractor is one of the stable 
critical points of the autonomous system. To study the stability of these critical points, we substitute 
linear perturbations x — ► x + Sx, y — > y + Sy, z — > z + Sz, u — ► u + Su and v — ► v + Sv about the critical 
point (x,y,z,u,v) into the dynamical system Eqs. (f3"5)) - (|4"2")) and linearize them. We get the evolution 
equations for the fluctuations as 



Sx' = (0i -3) Sx + xSOi (u5Q 2 + © 2 5u) + V2f 2 v 2 u 3 e^ 2 6Q 2 

-V2f 2 (3u 2 v 2 Su + 2u 3 vSv) Q^ 1 - 8C X , 



Sy' 


- (e 


Sz' 


= ( s 


Su' 


= u c 


Sv' 


= Sy, 



1 f2 -2 



-1 A 



5©i - V2xu@2 2 SQ 2 + V2 (xSu + uSx) 



01 + \f2xuQ- 1 ) Su 



(70) 
(71) 

(72) 

(73) 
(74) 



where 5©i, 502, SC±, SC 2l SQi and SQ 2 are the linear perturbations coming from ©i, ©2, C±, C 2l Q\ 
and Q 2 , respectively. The five eigenvalues of the coefficient matrix of the above equations determine the 
stability of the corresponding critical point. 

Now, we work out 5@i, 502, 5Ci, 5C2, 5Qi and 5Q2 one by one. From Eq. ([IS]) , we get 

(75) 
(76) 



5©i = 6exSx + 2ef 2 (u 2 vSv + uv 2 Su) + 4ry (y + v) (Sy + Sv) + 3jzSz. 
From Eq. (|45|) . we have 

502 = (A©2) 1 [zSz + exSx + ef 2 (u 2 vSv + uv 2 Su) +rj(y + v) (Sy + Sv)] 
From Eq. ([ST]]) , it is easy to find that 



V2a _ 
0C1 = zoz, 



SC 2 = —j= (xSz + zSx) . 
v2 



(77) 



Then, we obtain 5Qi and 5Q2 for Case (I) and Case (II) respectively, since they depend on the particular 
form of Q. For the Case (I), from Eq. (flTLj) . we obtain 



8 QI = - 3 MZ* l(Sy + Sv)-^L_Sz, 
2 V (y + vf yy ' V(y + v) 



For the Case (II), from Eq. it is easy to get 

SQi = -^(y + v)- 2 (5y + Sv) 
2r\ 



5Q 2 = --f3(z- 2 + l)Sz. 



SQ 2 = ~2 Crz 2 S Z - 



(78) 



(79) 
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We substitute the critical point (x, y, z, u, v) and ©i = as well as 82 given by Eq. (|58|) into Eqs. (|70|) - 
(j?4"|) with Eqs. (|75)) - l|75|) . The five eigenvalues of the coefficient matrix of these equations determine 
the stability of the corresponding critical point. We find that for both Case (I) Q = 2>(3Hp cw and 
Case (II) Q = 3aHp to t, the critical points can exist and are stable in proper parameter spaces, respec- 
tively. In other words, they are late time attractors. Needless to say, the particular parameter spaces for 
the existence and stability of these critical points are considerably involved and tedious. Since our main 
aim here is just to point out the fact that it can exist and is stable, we do not present those very involved 
expressions for the corresponding parameter space. Of course, one can work out with the help of Mathe- 
matica. Here, we just give several particular examples to support our statement. For the Example (1.1) 
of the Case (I) mentioned above, the corresponding eigenvalues are {—3.10805 + i 3.37335, —3.10805 — 
i 3.37335, -3.43776, -1.49442, -0.134622}; for the Example (1.2), the corresponding eigenvalues are 
{-3.10108 + ^3.12107, -3.10108 -i 3.12107, -2.95079, -1.67829, -0.108711}; for the Example (II.l) 
of the Case (II) mentioned above, the corresponding eigenvalues are {—2.68308 + i 2.6312, —2.68308 — 
i 2.6312, -2.42269, -1.6597, -0.032674}; for the Example (II.2), the corresponding eigenvalues are 
{-2.68735 + i 3.12848, -2.68735 - i 3.12848, -1.97246 + i 0.723607, -1.97246 -i 0.723607, -0.0785649}. 
Obviously, they are all stable. 



D. The first and second cosmological coincidence problems 

Here, following the argument in Sec. VI of [36|, we briefly show that the first and second cosmological 
coincidence problems can be alleviated at the same time in our model. As is well known, the approach 
most frequently used to alleviate the cosmological coincidence problem is the scaling attractor(s) in the 
dynamical system (see [34j, |36j, |60j, |6l|, [62J, [63J, [6J, |65|, [66J, 167J for examples). The most desirable feature of 
dynamical system is that the whole system will eventually evolve to its attractors, having nothing to do 
with the initial conditions. Therefore, fine-tuning is unnecessary. 

As is explicitly shown in this work, all stable attractors have the desirable properties, namely, their 
corresponding fi m and fl CVJ are comparable, while w CVJ < — 1. As mentioned in the end of Sec. IHII of the 
present paper, for the case of e and r\ have opposite signs, w cw can be larger than or smaller than —1. 
Thus, crossing the phantom divide is possible. For a fairly wide range of initial conditions with w CVJ > — 1, 
the universe will eventually evolve to the scaling attractor(s) with w cw < — 1, while the corresponding f2 TO 
and r2 cw are comparable. Thus, it is not strange that we are living in an epoch when the densities of dark 
energy and matter are comparable, and the EoS of dark energy is smaller than —1. In this sense, the first 
and second cosmological coincidence problems are alleviated at the same time in our model. 



V. CONCLUSION AND DISCUSSIONS 



In summary, the (generalized) Cheng- Weyl vector field coupled with a scalar field and its cosmological 
application are investigated in the present work. In our model, the dark energy is described as a mixture 
of a scalar field and a so-called "cosmic triad" of three mutually orthogonal Cheng- Weyl vector fields. We 
derive the effective energy density and pressure of the "mixed" dark energy, and the equations of motion 
for the scalar field and the Cheng- Weyl vector field, respectively, by using the spatial volume-average 
procedure. The cosmological evolution of this "mixed" dark energy is studied. We find that the effective 
EoS can cross the phantom divide Wd e = — 1 in the case of e and n having opposite signs. The first and 
second cosmological coincidence problems can be alleviated at the same time in our model. On the other 
hand, it is easy to see that all stable attractors have w e f / — — 1 . Although the EoS of dark energy u> cw 
can be smaller than —1, the big rip never appears in this model. The fate of our universe is an inflationary 
phase, in which the Hubble parameter H = H = const. 

we finish this paper with some remarks. In this work, to cross the phantom divide, e and 77 should have 
opposite signs. This means that either the scalar field or the vector field has a negative sign of its kinetic 
term. As is argued in [7l|, these ghost-type fields are possible and viable. Second, one might notice that 
from Eqs. (|6"3"|) . l|o3|) . (|6"5f and (|6"9"|) . at the attractors, the appropriate value VL m ~ 0.25 at [3 ~ 1/3 or 
a ~ 1/4 results w cw ~ —4/3, for 7 = 1. However, this is not inconsistent with the observations [72l |. 
The situation becomes more comfortable when one is aware of the possibility that the universe perhaps 
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nears but has not reached the attractors. Finally, at the attractors, fl m is constant and not equal to zero, 
thanks to the interaction between the vector fields and the background matter. In fact, the interaction 
between dark energy and dark matter can be constrained by the cosmological observations, see [37], [ill for 
examples. However, there is still a long way to obtain some strict constraints on this interaction. 
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